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ABSTRACT 

Period-colour (PC) and amplitude-colour (AC) relations are studied for the Large 
Magellanic Cloud (LMC) Cepheids under the theoretical framework of the hydrogen 
ionization front (HIF) - photosphere interaction. LMC models are constructed with 
pulsation codes that include turbulent convection, and the properties of these models 
are studied at maximum, mean and minimum light. As with Galactic models, at max- 
imum light the photosphere is located next to the HIF for the LMC models. However 
very different behavior is found at minimum light. The long period {P > lOdays) LMC 
models imply that the photosphere is disengaged from the HIF at minimum light, sim- 
ilar to the Galactic models, but there are some indications that the photosphere is 
located near the HIF for the short period {P < 10 days) LMC models. We also use the 
updated LMC data to derive empirical PC and AC relations at these phases. Our nu- 
merical models are broadly consistent with our theory and the observed data, though 
we discuss some caveats in the paper. We apply the idea of the HIF-photosphere in- 
teraction to explain recent suggestions that the LMC period-luminosity (PL) and PC 
relations are non-linear with a break at a period close to 10 days. Our empirical LMC 
PC and PL relations are also found to be non-linear with the F-test. Our explanation 
relies on the properties of the Saha ionization equation, the HIF-photosphere interac- 
tion and the way this interaction changes with the phase of pulsation and metallicity 
to produce the observed changes in the LMC PC and PL relations. 
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1 INTRODUCTION 

ICodd (^ij) found that the Galactic Cepheids follow a spec- 
tral type that is independent of their pulsational periods at 
maximum lig ht and gets later as the periods increase at min- 
imum light. iSimon et all il993l hereafter SKM) used radia- 
tive hydrodynamical models to explain these observational 
phenomena as being due to the location of the hydrogen 
ionization front (HIF) relative to the photosphere. Their re- 
sults agreed very well with Code's observation. SKM further 
used the Stefan-Boltzmann law applied at the maximum and 
minimum light, togeth er with th e fact that radial variation 
is small in the optical JCoxlll980ri . to derive: 
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log Trriax log T^min — (Vmin Vmax), (1) 

where Tyy^ax/min are the effective temperature at the max- 
imum/minimum light, respectively. If Tmax is independent 
of the pulsation period P (in days), then equation (1) pre- 
dicts there is a relation between the l/-band amplitude and 
the temperature (or the colour) at minimum light, and vice 
versa. In other words, if the period-colour (PC) relation 
at maximum (or minimum) light is flat, then there is an 
amplitude-colour (AC) relation at minimum (or maximum) 
light. Equation (1) has shown to be valid theoretically and 
obser vationally for the classica l Cepheids and RR Lyrae va ri- 
ables jKanbur fc Ngeowll2004 iKanbur fc Ferrianddl2005D . 

For the RR Lyrae variBbles~ lKanbuil i^Qsli and 
iKanbur fc Phillipj il996l) used linear and non-linear hy- 
drodynamic models of RRab stars in the Galaxy to ex- 
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plain why RRab s tars follow a flat PC relati on at min- 
imum light. Later, iKanbur fc Fernando! ^2005^ used MA- 
CHO RRab stars in the LMC to prove that LMC RRab 
stars follow a relation such that higher amplitude stars are 
driven to cooler temperatures at maximum light. Similar 
studie s were also carried out fo r Cepheid variables , as in 
SKM, iKanbur fc PhiU pj lHooil. iKanbur^fc Ngeowl J2004 
hereafter Paper I) and IKanbur et al J ('2004', hereafter Paper 
II). In contrast to the RR Lyrae variables, Cepheids show a 
flat PC relation at the maximum light, and there is a AC re- 
lation at the minimum light. Therefore, the PC relation and 
the AC relation are intimately connected. All these studies 
are in accord with the predictions of equation (1). 

In Paper I, the Galactic, Large Magellanic Cloud (LMC) 
and Small Magellanic Cloud (SMC) Cepheids were analyzed 
in terms of the PC and AC relations at the phase of maxi- 
mum, mean and minimum light. One of the motivations for 
this paper originates from recent studies on the non-linear 
LMC PC relation (as well as the period-luminosity, PL, re- 
lation. See Paper I: to.mmann et al.l 120021 : ISandage et alJ 
l2004l:lNgeow et al.ll2005t) : the optical data are more consis- 
tent with two lines of differing slopes which are continuous 
or almost continuous at a peri od close to 10 days. Paper 
I also applied the the F-test iWeisberd figSOl) to the PC 
and AC relations at maximum, mean and minimum V-band 
light for the Galactic, LMC and SMC Cepheids. The F-test 
results implied that the LMC PC relations are broken or 
non-linear, in the sense described above, across a period of 
10 days, at mean and minimum light, but only marginally 
so at maximum light. The results for the Galactic and SMC 
Cepheids are similar, in a sense that at mean and minimum 
light the PC relations do not show any non-linearity and 
the PC (max) relation exhibited marginal evidence of non- 
linearity. For the AC relation, Cepheids in all three galaxies 
supported the existence of two AC relations at maximum, 
mean and minimum light. In addition, the Cepheids in these 
three galaxies also exhibited evidence of the PC- AC connec- 
tion, as implied by equation (1), which give further evidence 
of the HIF-photosphere interactions as outlined in SKM. 

To further investigate the connection between equation 
(1) and the HIF-photosphere interaction, and also to explain 
Code's observations with modern stellar pulsation codes. 
Galactic Cepheid models were constructed in Paper II. In 
contrast to SKM's purely radiative models, the stellar pul- 
sation codes used in Paper II in cluded the treatmen t of tur- 
bulent convection as outlined in lYecko et alJ lll998f) . One of 
the results from Paper II was that the general forms of the 
theoretical PC and AC relation matched the observed rela- 
tions well. The properties of the PC and AC relations for 
the Galactic Cepheids with log(P) > 0.8 can be explained 
with the HIF-photosphere interaction. This interaction, to 
a large extent, is independent of the pulsation codes used, 
the adopted ML relations, and the detailed input physics. 

The aim of this paper is to extend the investigation 
of the connections between PC- AC relations and the HIF- 
photosphere interactions in theoretical pulsation models of 
LMC Cepheids, in addition to the Galactic models presented 
in Paper II. In Section 2, we describe the basic physics of 
the HIF-photosphere interaction. The updated observational 
data, after applying various selection criteria, that used in 
this paper are described in Section 3. In Section 4, the new 
empirical PC and AC relations based on the data used are 



presented. In Section 5, we outline our methods and model 
calculations, and the results are presented in Section 6. Ex- 
amples of the HIF-photosphere interaction in astrophysical 
applications are given in Section 7. Our conclusions fc dis- 
cussion are presented in Section 8. Throughout the paper, 
short and long period Cepheid are referred to Cepheids with 
period less and greater than 10 days, respectively. 



2 THE PHYSICS OF HIF-PHOTOSPHERE 
INTERACTIONS 

The partial hydrogen ionization zone (or the HIF) moves 
in and out in the mass distribution as the star pulsates. It 
is possible that the HIF will interact with the photosphere, 
defined at optical depth (r) of 2/3, at certain phases of pul- 
sation. For example, SKM suggested that this happened at 
maximum light for the Galactic Cepheids, as the HIF is so 
far out in the mass distribution that the photosphere occurs 
right at the base of the HIF. The sharp rise of the opac- 
ity wall (where the mean free path goes to zero) due to the 
existence of HIF prevents the photosphere moving further 
into the mass distribution and hence erases any "memory" 
of global stellar conditions, including the underlying PC re- 
lation. This lead to a flat relation between period & temper- 
ature, period fc colour and period & spectral type at maxi- 
mum light, as seen in SKM and Paper II. At other phases, 
since the HIF does not interact with the photosphere, the 
temperature of the star (or the colour) follows the underly- 
ing global PC relation. 

The HIF-photosphere interaction also relies on the 
properties of the Saha ionization equation and the struc- 
tural properties of the outer envelopes of Cepheids. It is 
well known that the partition functions in the Saha ion- 
ization equation are formally divergent unless some atomic 
physics is used to truncate them. In the pulsation codes 
we used, we approximate the partition functions of various 
atoms by their ground state statistical weights. The proper- 
ties of the Saha ionization equation in Cepheid envelopes are 
such that hydrogen starts to ionize at a temperature that 
is almost independent of density, for a certain range of low 
densities. Outside of this range of density, the density depen- 
dence increases. Thus, when the photosphere is very close to, 
or engaged with the HIF and the density of these regions is 
reasonably low, the temperature of the photosphere is less 
dependent on the surrounding density and hence the global 
stellar parameters. At higher densities, the temperature at 
which hydrogen ionizes becomes more sensitive to density 
and hence more sensitive to global stellar parameters. 

If the photosphere is far from the HIF, or disengaged, 
then the location of the photosphere and hence the tem- 
perature of the photosphere, is again strongly dependent on 
density and hence on global stellar parameters. That is why 
the photosphere needs to be close to, or engaged with the 
HIF for this effect to take place. Moreover, this dependence 
on density is not sharp so that for " low" and " high" densities 
the density dependence of the photospheric temperature is 
w eak and strong respe ctively. An examination of figure 15.1 
in ICox fc Giuli lll968h demonstrates that this is plausible. 
Thus as the star pulsates, the photospheric temperature has 
a density dependence that can be strong or weak depending 
on phase. An example where the density dependence is weak 
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are the Galactic long period Cepheids at maximum light 
(SKM, Paper II): these Cepheids display a flat PC relation 
at maximum light. These properties of the HIF-photosphere 
interaction can, in turn, affect the temperature of the pho- 
tosphere and hence the colour of the Cepheid. 

Here we investigate the idea that LMC Cepheids with 
periods below 10 days are such that the HIF and photo- 
sphere are engaged through most of the pulsation cycle. At 
periods greater than 10 days, the photosphere only engages 
with the HIF at maximum light. The transition is sharp be- 
cause the photosphere is either at the base of the HIF or 
it is not. The transition occurs because as the period in- 
creases, the L/M ratio increases and this implies the HIF is 
located further inside in the mass distribution, chan ging the 
phase at which it can interact with the photosphere (|KanbuJ 
Il995r . The structure of Galactic Cepheids is such that this 
interaction only occurs at maximum light, even for Cepheids 
with periods shorter than 10 days. 
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3 UPDATED LMC CEPHEID DATA 

In Paper I, we constructed the light curves of fundamental 
mode Cepheids in the LMC by using the extensive photo- 
metric dataset in the OGLE (Optical Gravitational Lensing 
Experiment) database. However, the dataset used in Paper I 
was downloaded in 2002, prior to the updated version of the 
dataset that was available after April 24, 2003 (OGLE web- 
site, Udalski 2004 [private communication]). The updated 
version includes additional V- and /-band data for most 
of the Cepheids. In addition, the periods have been refined 
by the OGLE team using the complete set of photomet- 
ric data. Due to these r easons, we decided to repeat the 
light curve construction iNgeow et aljl20o3l with the up- 
dated data and periods. Since the Cepheids in the OGLE 
database are truncated at log(P) ~ 1.5, due to the sat- 
uration of the CC D detector for the lon ger period (hence 
brighter) Cepheids llUdalski et alJll999a|). we include sorn e 
additional LMC Ce pheid data f rom 'Mo ffett et alJ Jl99d) . 
iBarnes et al] lll999l) and ISebo e t al. ( 200l) to extend the 
period coverage to log(P) > 1.5 in our sample. The require- 
ments that govern our choice of the published photometric 
data are: (a) latest observations that use the modern day 
CCD cameras; (b) high quality data with large number of 
data points per light curve, which provide uniform phase 
coverage and small scatter of the light curve; and (c) as 
homogeneous as possible (i.e., from a minimal number of 
sources) to avoid any additional systematic errors. These 
requirements are essential to construct accurate light curves 
to allow the estimation of colours and magnitudes at maxi- 
mum, mean and minimum light for our PC and AC study. 
Hence we did not include some of the older photometric data 
in this study. 

The photometric data of all Cepheids, com pris- 
ing 7 7 1 from OGL E database, 14 from iMoffett et alJ 
1 1998 1) -t jBarnes et al.l lll999l) and 39 from ISebo et alJ 
I 2002 ) ^ , were mainly fit with n = 4 to n = 8 Fourier expan- 



Figure 1. The PC relation at maximum (upper panel) and mini- 
mum (lower panel) light for the long period Cepheids in our sam- 
ple (after removal of outliers). The solid and dashed lined arc the 
regressions that exclude and include the longest period Cepheids. 
It is clear that including the Cepheids with log(P) > 1.8 has made 
the PC(max) relation becomes steeper, even though the two re- 
gressions are consistent to each others at the minimum light. 



sions (n is the order of Fourier expansion) using the sim- 
ulated annealing method described in Ngcow ct al ( 2003) 
to the V- and /-band photometric data. This is in contrast 
to Paper I that only applied n = 4 Fourier fits. However, 
for some of the OGLE long period Cepheids (P > 11.5 
days), it was found out that the quality of the fitted light 
curves could be improved by using a higher order Fourier 
expansion, hence we extended the fit to n = 12 for these 
long period Cepheids. All the fitted light curves were visu- 
ally inspected and the best-fit light curves from the differ- 
ent orders of the Fourier expansions were selected. To the 
best of our knowledge, this analysis also represents a ma- 
jor improvement in the Fourier analysis of the OGLE data. 
The extinction is corrected with the standard procedure, 
i.e. {V - I)o = (V - /) - [Ry - Ri) X E(B - V) with 
Rv = 3.24 and Ri = 1.96 iUdalski et alJll999bl) . The val- 
ues of E{B — V) for each OGLE Cepheids are taken from 
the OGLE database JUdalski et all Il999bft. while fo r the 
Cepheids in 'Moffett et al.1 Jl998t) -l jBarnes et aP and 
in Scbo ct al. (2002), the values o f E(B — V) are ado pted 
from lSandage et al.l ll20o3) and/or lPersson et alJ ^Qoi). 

To guard against some "bad" Cepheids or other contam- 
ination in our sample, and select only the good Cepheids 
in both bands, we removed some Cepheids in the sample 
accor ding to the following criteria (see also iKanbur et aP 
l2003h: 



^ For ISebo et alj l2002^ data, the photometry data are only 
available for the long period Cepheids. To avoid duplic ation, 
we remove the long period Cepheids i n ISe bp et al.l J2002ft that 
are labeled "Ogle", those ISebo et al.l l20o3r ^dentified as the 



(i) Cepheids without V- and/or /-band photometry, or 
the number of data per light curve (in either bands or both) 
is too low to fit a 71 = 4 Fourier expansion. 



OGLE Cepheids, a nd 5 "HV" Cepheids that are t he same in 
iMoffett et ai]^199g^ . This left 39 "HV" Cepheids from lSebo et alJ 
i2002l) sample. 
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(ii) Cepheids with poorly fitted or unacceptable V- 
and/or /-band light curves in the sample, such as those with 
a large scatter of data points or with bad-phase coverage 
(large gaps between the phased data points). Most of the 
magnitudes, as well as the colours, at the maximum and/or 
minimum light from these fitted light curves are very uncer- 
tain. 

(iii) Cepheids with possible duplicity in the OGLE sam- 
ple. Some of the possible duplicated Cepheids were removed 
in the O GLE database by consulting table 4 of lUdalski et alJ 
Jl999d) . 

(iv) Cepheids with unusual colour. We first plot out (as 
in Figure |2Ia]) the extinction corrected PC relation at mean 
light. The plot shows that there are number of outliers in the 
period-colour plane, mostly with log(P) < 1.0. The presence 
of these outliers is probably due to: (a) their extinction is 
either over- or under-estimated; (b) they have blue or red 
companions that cannot be resolved due to the problems 
of blending; or (c) other unknown physical reasons. A de- 
tailed investigation of these outliers is beyond the scope of 
this paper, but it is clear that they should be removed from 
the sample. These outliers are removed with the adopted 
colour-cut of 0.35 < {V - 7)5''='"' < 0.95, a compromise be- 
tween maximizing the number of Cepheids in the sample 
and excluding the Cepheids with unusual colour^. 

(v) Cepheids with unusually low (or high) amplitude. 
Some Cepheids with unusually low V- and /-band ampli- 
tudes were found in the sample. Their amplitudes are typ- 
ically 2^3 times smaller as compared to the amplitudes 
of other Cepheids at given period. Some examples of the 
light curves for thes e low amplitude Cepheids are given in 
iKanbur et alJ ll2003t) . In addition, most of the light curves 
for these low amplitude Cepheids can be fitted with n = 4 
Fourier expansion, while other Cephei ds with "normal" am - 
plitude may require higher order fits. iKanbur et al.l ll2003ll 
has briefly discussed some possible physical reasons for these 
Cepheids to have such low amplitudes, e.g. they are just 
entering or leaving the fu ndamental mode instability strip 
jBuchler fc KoUathl liooi) or they have different chemical 
composition (see, e.g., [Paczvriski & Pindor 200^). The de- 
tailed investigation of these low amplitudes Cepheids is be- 
yond the scope of this paper. Here, we apply a conserva- 
tive amplitude cut of 0.3 mag. in the V-band to remove 
the low amplitude Cepheids. Besides that, we also remove 
OGLE-286532 (with unusually low amplitude) and HV-2883 
(with unusually high amplitude) as they are c lear outliers i n 
the log(P )-amplitude pl ot (not shown, but see lNgeowl2005l) . 
Note that iPietrzviiski et al. (2004|) applied a cut of 0.4 mag. 
to remove the low amplitude Cepheids in NGC 6822. Other 
examples o f removing the low-am plitude Cepheids can also 
be found in lWavman et all 1^8^. 

(vi) Cepheids with log(P) < 0.4 and log(P) > 1.8. 
In order to guard against possible contaminati on from 
the first overtone Cepheids lUdalski_.e,t al."l999a') and to 
be consistent with the previous studies IlKanbur fc Ngeowl 

^ Note that this colour-cut is only applied to the short period 
Cepheids because there are not many long period Cepheids in 
the sample. Furthermore, the PC (mean) relation for the long 
period Cepheids is steeper than the short period Cepheids, hence 
as the period increases the colour becomes redder, which can be 
excluded from the colour-cut. 



l2004l:ISanda'ge' et al .120041 : lUdalskT et al .ligggaTi . we removed 
Cepheids with log (P) < 0.4 (see further justification in 
iNgeow et af]|2005^ . Regarding the removal of Cepheids with 
log(P) > 1.8, a preliminary analysis of the PC relation re- 
veals that few of the longest Cepheids should be removed 
from the sample, because they are clear outliers in the PC 
plot at maximum light (see upper panel of Figure 0. With- 
out these longest period Cepheids, the PC (max) relation for 
the long period Cepheids is fiat, which is consistent with 
the results found in Paper I. The hypothesis of the HIF- 
photosphere interaction also suggests the fiatness of the PC 
relation at maximum light for long period Cepheids. How- 
ever, as the period gets longer (with log[P1 > 1.5), th e 
photosphere disengages from the HIF llSimon et alj[l993) . 
These longest period Cepheids have biased the slope of the 
PC(max) relation by making the slope becomes steeper. 

These selection criteria are guided mainly by the philoso- 
phy that it is better to lose some "bad" but real Cepheids 
rather than inc luding those sp urious and doubtful Cepheids 
in the sample llLeonard et alj. 20031 . or those with bad fit- 
ted light curves that will give inaccurate measurements of 
the maximum and minimum light. Hence, the final sam- 
ple consists of 641 LMC Cepheids that will be considered 
further. The locations of the outliers from various selec- 
tion criteria are shown in Figure |5| for the PC(mean) re- 
lation, V'-band PL relation, P2i-log(P) relation'' and the 
colour- magnitude diagram (CMD). Note that some of the 
outliers are located within the "good" Cepheids. However 
they can be eliminated due to various physical reasons as 
given above, especially those with poorly fit light curves 
that will give inaccurate measurements at maximum, mean 
and minimum light. A simple sigma-clipping algorithm (e.g., 
iUdals ki ct al. 1999a) wil l not be able to remove these out- 
liers iKanbur et al.li2003l) . 



4 THE NEW EMPIRICAL PC AND AC 
RELATIONS 

To construct the empirical PC fc AC relations, we used the 
following quantities from the Fourier fits to the Cepheid data 
as obtained from previous section: 

• l/-band amplitude: the difference of the numerical max- 
imum and minimum from the Fourier expansion, Vamp = 

• (V - /)"""^: defined as Vmax - Iphmax, where Iphmax is 
the /-band magnitude at the same phase as Vmax- 

• (V - /)™=<'": defined as Ao{V) - Ao{I), w here Ap is the 
mean value from the Fourier expansion (see iNgeow et alJ 
l2003h . This is very similar to the conventional definition of 
the mean colour, < V > — < I > , where <> denotes the 
intensity mean. 

• {V - 7)P'i™<=''": defined as V^ean - Iphmean, where 
Iphmean is the /-band magnitude at the same phase as 
Vmean- Vjnean is the V-baud magnitude closest to A(i{V), 
the mean value from Fourier expansion. 

• {V - /)™'": defined as Vmin - /phmin, where Iphmin is 
the /-band magnitude at the same phase as Vmin- 

^ R21 = A2/A-1, wher e An are the Four ier amplitudes. See 
ISimon fc Le^ il98ll) and lNgeow et al.l ^2003^ for details. 
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Figure 2. -The various relations and the colour-magnitude diagram (CMD) for the LMC Cepheids, after removing Cepheids with 
log(P) < 0.4 and those lacking of V-band photometry. The symbols are (see text for details): filled circles = Cepheid with poor fit of 
the light curves; open squares = the rejected Cepheids after colour-cut of 0.35 < {V — /)™°^" < 0.95; open triangles = the excluded 
Cepheids after the amplitude cut; filled triangles = the removed Cepheids for possible duplicity; and open circles = the longest period 
Cepheids, which are the outliers in PC (max) relation. The dots are the remaining good Cepheids in our sample. 



These quantities have been corrected for extinction as men- 
tioned in previous section. The empirical LMC PC and AC 
relations at maximum, mean and minimum light for all, long 
and short period Cepheids are summarized in Table 
and the corresponding plots are presented in Figure 01 
respectively. 

To test the non-linearity of the PC and AC relations, 
or the "break" at a period of 1 days, we apply t he F- 
test as given in Paper I and in iNgeow et al.l (|2TO^ . The 
null hypothesis in the f-test is single line regression is suf- 
ficient, while the alternate hypothesis is that two lines re- 
gressions with a discontinuity (a break) at 10 days is nec- 
essary to fit the data. The probability p{F), under the 
null hypothesis, can be obtained with the corresponding 
_F-values and the degrees of freedom. In general, the large 
value of _F (equivalent to the small value of p[F]) indi- 



cates that the null hypothesis can be rejected. For our 
sample, F ~ 3.0 when p{F) = 0.05 (the 95% confident 
level), therefore the null hypothesis can be rejected if the 
_F-value is greater than 3 with more than 95% confident 
level and the data is more consistent with the two-line re- 
gression. A glance of Tableland Figure |3| suggests that the 
LMC PC relations are broken at maximum, mean and mini- 
mum light. These are confirmed with the F-test results with 
-Fpc (max, mean, phmean, min) — {9.49, 8.59, 16.8, 11.8}. 
Similarly, the _F-test results for the AC relation are: 
^AC (max, mean, phmean, min) = {139, 162, 157, 171}. 
Hence, the LMC PC and AC relations are non-linear (hence 
broken) at maximum, means and minimum light. Note that 
the flatness of the long period PC(max) relation as given in 
Table0 (0.012 ± 0.069) is in good agreement with the slope 
found in Paper I (—0.031 ± 0.101). Recall that equation (1) 
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Table 1. The LMC period-colour relation in the form of (V — /) = 
alog(P) + 6, and cr is the dispersion of the relation. 



Phase 


a 




b 




f7 




A u ^^ — 










Maximum 


0.168 ±0.017 


0.328 


± 


0.012 


0.099 


Mean 


0.228 ±0.013 


0.479 


± 


0.009 


0.075 


Phmean 


0.248 ±0.014 


0.489 


± 


0.010 


0.081 


Minimum 


0.297 ±0.013 


0.532 


± 


0.009 


0.075 




Long, N 


= 63 








Maximum 


0.012 ±0.069 


0.502 


± 


0.086 


0.098 


Mean 


0.284 ±0.053 


0.428 


± 


0.066 


0.075 


Phmean 


0.383 ± 0.065 


0.347 


± 


0.081 


0.092 


Minimum 


0.311 ±0.056 


0.539 


± 


0.070 


0.081 




Short, N 


= 578 








Maximum 


0.263 ±0.029 


0.272 


± 


0.018 


0.097 


Mean 


0.153 ±0.022 


0.523 


± 


0.014 


0.074 


Phmean 


0.141 ±0.023 


0.552 


± 


0.014 


0.078 


Minimum 


0.210 ±0.022 


0.582 


± 


0.014 


0.073 



Table 2. The LMC amplitude-colour relation in the form of {V — 
I) = aVamp + b, and cr is the dispersion of the relation. 



Phase 


a 




b 




a 




All, N = 


641 








Maximum 


-0.274 ±0.019 


0.641 


± 


0.015 


0.093 


Mean 


0.029 ±0.019 


0.611 


± 


0.014 


0.092 


Phmean 


0.074 ± 0.021 


0.602 


± 


0.016 


0.099 


Minimum 


0.159 ±0.020 


0.616 


± 


0.015 


0.098 




Long, N 


= 63 








Maximum 


-0.209 ± 0.049 


0.716 


± 


0.048 


0.086 


Mean 


0.148 ±0.048 


0.637 


± 


0.047 


0.085 


Phmean 


0.230 ±0.058 


0.600 


± 


0.057 


0.103 


Minimum 


0.240 ±0.046 


0.693 


± 


0.045 


0.082 




Short, N -- 


= 578 








Majcimum 


-0.419 ±0.019 


0.729 


± 


0.014 


0.076 


Mean 


-0.130 ±0.018 


0.708 


± 


0.013 


0.074 


Phmean 


-0.099 ±0.109 


0.708 


± 


0.014 


0.078 


Minimum 


-0.008 ± 0.020 


0.716 


± 


0.014 


0.079 



predicts that if the PC relation is flat at maximum hght, 
then there is a correlation between the amplitude and the 
colour at minimum lieht. This is seen in Table (and m 
Figure O for the long period AC(min) relation, with a slope 
of 0.240 ± 0.046. 



5 METHODS AND LMC MODELS 

The stellar pulsa tion codes we u sed are both linea r 
jYecko et alJ Il998t) and non-linear feoUath et al.l l2002l) . 

These codes, which inc lude a 1-D turbulent convection 
recipe jYecko et al.lll99 8<l. are the same as in Paper II. 
Briefly speaking, the codes take the mass (M), luminosity 
(L), effective temperature {T^ff) and chemical composition 
{X, Z) as input parameters. The chemical composition is 
set to be {X,Z) = (0.70,0.008) to represent the LMC hy- 
drogen and metallicity abundance (by mass) . The mass and 
luminosity are obtained from the ML relations calculated 
from evolutionary models. The T^f / are chosen to ensure the 
models oscillate in the fundamental mode and located inside 
the Cepheid instability strip. The pulsation periods for the 



models are obtain ed from a linear non-adiabatic analysis 
jYecko et alJll998^ . All other parameters used in the pulsa- 
tion codes had the same values for the LMC and Galactic 
models (Paper II). This included the a parameters that are 
part of the turbulent convection recipe, though see Section 
8. Of course, one variable parameter was the metallicity. 
The only other difference between this study and Paper II, 
besides the metallicity, is the value set for the artificial vis- 
cosity parameter, Cq. In this study, we set Cq — 16.0 for the 
LMC models to improve the shape of the theoretical light 
curves, in contrast to the value of 4.0 used for the Galactic 
models. 

In Pa per II, the ML rela tions are adopted from IChiosil 
ljl989h and lBono et~ai] i200(tl . In order to be consistent with 
previous work, the ML relations used in this paper w i ll also 
be adopted from these two sources. However. IChiosil (^8^ 
only provided two ML relations, one for Z = 0.020 which are 
used in Paper II, and another one for Z = 0.001. Hence we 
have to adopt the second ML relation for the LMC models. 
Even though the LMC metallicity is higher than Z = 0.001, 
the LMC is still conside red as a low metallicity system in 
the literature. Hence the IChiosi il989) ML relation can be 
approximately applied for the LMC models. An anonymous 
referee pointed out that an interpolation of the lChiosil (Il989l) 
ML relations between Z = 0.020 and Z = 0.001 should 
also be used. We have included the interpolated ML rela- 
tion in our model calculations. In the context of the HIF- 
photosphere interaction, it is the ML relation which dictates 
at what period and at what phases this will occur. Stellar 
evolutionary theory changes the ML relation as a function 
of metallicity. Hence the coefflcients of the ML relation are 
important in determining the nature of the HIF-photosphere 
interaction (Paper II). In short, the ML relations used are: 

(i) ML relation given in lBono etall J2000l) : 

log(L) = 0.90 + 3.351og(M) -I- 1.361og(y) - 0.341og(Z), 
= 3.35 log(M) + 0.886. (2) 

(ii) ML relation given in lChiosi Jl989l) : 

log(L) = 3.22 log(M) + 1.511. (3) 

(iii) ML relation interpolated between two IChiosil (^8^ 
relations at Z = 0.02 and Z = 0.001 to yield a relation at 
Z = 0.008: 

log{L) = 3.36 log(Af) -1- 1.295. (4) 

The units for both A4 and L are in Solar units. Note that 
these ML relations cover reasonably broad L/AI ratios given 
in the literature. The input parameters for the LMC models 
with these ML relations and the periods calculated from 
linear non-adiabatic analysis are given in Table |H] 

After the full amplitude models are constructed from 
the pulsation codes, the temperature and the opacity profile 
can be plotted in terms of the internal mass distribution 
(log[l — Mr/M], where Mr is mass within radius r and M 
is the total mass) at a given phase of pulsation. As in Paper 
II, the locations of the HIF (sharp rise in the temperature 
profile) and photosphere (at optical depth r = 2/3) can be 
identified in the temperature profile. To quantify the HIF- 
photosphere interaction (if the photosphere is next to the 
base of the HIF or not, see also Paper II), we calculate the 
"distance", A, in log(l - Mr/M) between the HIF and the 
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Figure 3. The period-colour (PC) relations for the LMC Cepheids at maximum, means and minimum light. The open and filled circles 
are for short and long period Cepheids, respectively. The solid and dashed lines are the fitted PC relations for the short and long period 
Cepheids, respectively. 




Figure 4. Same as Figure |3 but for the amplitude-colour (AC) relations. 



photosphere from the temperature profile. The definition of 
A can be found in Paper II. A small A means there is a 
HIF-photosphere interaction, and vice versa. 

The theoretical quantities from the models can be com- 
pared to the observed quantities using the following pre- 
scriptions: 

(i) As in Paper II. we use the BaS eL atmosphere database 
jLeieunell2002l:IWestera et alJl2002ll to construct a fit giving 
temperature and effective gravity as a function of {V — I) 
colour. The effective gravity is obtained at the appropriate 
phase from the models (see Paper II). These prescriptions 
are used to convert the temperatures to the {V — I) colours. 
The bolometric corrections (BC) are obtained in a similar 
manner. The anonymous referee has suggested that {V — I) 



may not be a good way to convert between temperature and 
colour unless both of the micro-turbulence and surface grav- 
ity are included. As indicated above this is the case, and 
in any case our results and those of Paper II for Galactic 
models, show good agreement between theory and observa- 
tions. A number of previous authors have used this method 
and some authors commented that th is colour can be used 
as an indicator of te mperature (e.g. iBeaulieu et al.lf200ll : 
[Pammann et aT] |2003l) . The empirical relations we studied 
in this series were also mainly in the {V — I) colour. 

(ii) In addition to t he BaSeL atmospher e, we also use 
the atmosphere fit from lSandage et alJ (|l99tt) . referring this 
as the SBT atmosphere in our paper. The SET atmosphere 
does include both of the effective gravity and the micro- 
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Table 3. Input parameters for LMC Cepheid models with periods 
obtained from a linear analysis. The periods, Po and Pi, are re- 
ferred to the fundamental and first overtone periods, respectively. 
Similarly for the growth rate, rj. Both of the mass and luminosity 
are in Solar units, the temperature is in K and the period is in 
days. 



Table 4. Temperatures at maximum and minimum light from 
full-amplitude non-linear model calculations. The periods, lumi- 
nosity and temperature are in days, Lq and K, respectively. 



M 


log(L) 




Po 


VO 


Pi 


m 




ML Relation from Bono et al. 


C2000') 




11.0 


4.375 


5050 


46.4155 


0.124 


28.98 


-0.118 


10.0 


4.236 


5100 


35.6727 


0.091 


22.92 


-0.093 


9.50 


4.161 


5250 


28.2406 


0.094 


18.68 


-0.046 


9.10 


4.099 


5260 


25.3960 


0.082 


16.92 


-0.042 


8.75 


4.042 


5310 


22.3804 


0.076 


15.07 


-0.027 


8 40 


3.982 


5380 


19.3886 


0.071 


13.20 


-0 008 


7.95 


3.902 


5330 


17.7750 


0.055 


12.09 


-0.027 


7.00 


3.717 


5410 


12.6085 


0.035 


8.722 


-0.018 


6.55 


3.620 


5490 


10.2940 


0.031 


7.183 


-0.006 


6.40 


3.587 


5485 


9.81474 


0.027 


6.853 


-0.010 


6.00 


3.493 


5510 


8.37226 


0.020 


5.866 


-0.014 


5.90 


3.468 


5500 


8.12498 


0.017 


5.691 


-0.018 


5.80 


3.443 


5525 


7.69466 


0.017 


5.400 


-0.015 


5 70 


3.418 


5560 


7.23505 


0.017 


5.090 


-0 009 


5.30 


3.312 


5600 


6.01283 


0.012 


4.244 


-0.009 




ML Relation from C 


Chiosi C1989') 




7.20 


4.272 


5380 


40.2561 


0.275 


24.51 


-0.162 


6.80 


4.192 


5380 


35.4374 


0.264 


21.91 


-0.122 


6.20 


4.063 


5410 


28.2629 


0.225 


17.94 


-0.076 


5.95 


4.005 


5420 


25.6378 


0.211 


16.43 


-0.060 


5.40 


3.869 


5510 


19.4314 


0.170 


12.82 


-0.015 


5.15 


3.803 


5510 


17.5523 


0.160 


11.66 


-0.007 


4.65 


3.660 


5490 


14.3143 


0.131 


9.611 


-0.010 


4.20 


3.518 


5510 


11.3659 


0.101 


7.729 


-0.011 


4.00 


3.450 


5545 


10.0011 


0.089 


6.854 


-0.005 


3.95 


3.432 


5540 


9.77393 


0.085 


6.701 


-0.008 


3.80 


3.378 


5550 


8.94637 


0.075 


6.157 


-0.009 


3.70 


3.341 


5575 


8.31297 


0.070 


5.745 


-0.005 


3.65 


3.322 


5570 


8.10751 


0.066 


5.605 


-0.008 


3.60 


3.302 


5530 


8.09994 


0.058 


5.583 


-0.022 


3.60 


3.302 


5600 


7.71463 


0.065 


5.352 


-0.001 




ML Relation from interpolated Chiosi ^^89) 




6.80 


4.092 


5280 


30.9661 


0.207 


19.51 


-0.091 


5.20 


3.701 


5340 


15.9744 


0.094 


10.63 


-0.050 


4.40 


3.457 


5460 


10.0562 


0.056 


6.898 


-0.028 


4.20 


3.389 


5550 


8.51317 


0.054 


5.906 


-0.007 


3.80 


3.243 


5630 


6.46303 


0.039 


4.532 


-0.002 



turbulence in their table 6 for the temperature and colour 
conversion. These conversions are tabulated for two micro- 
turbulence velocities of 1.7kms~^ and 5.0kms~^, as well as 
for various metallicities. To apply these conversions to our 
LMC models, we first interpolated the conversions between 
[A/H] = 0.0 and [A/H] = -0.5 to [A/H] = -0.3, which is 
appropriate for the LMC metallicity. The {V — I) colours at 
the maximum, mean and minimum Hght are then obtained 
from the given effective temperature and the effective grav- 
ity for both of the micro-turbulence velociti es. 

(iii) We use the prescriptions given in iBeaulieu et alJ 
l|200ll) to convert the observed colours to the temperatures 
appropriate for the LMC data as follows: 

log(3) = 2.62-L211og(P), 
log(re//) = 3.91545 -f 0.0056 log(g)- 0.2487(1/ -/)o, 
AT = log(r,//) - 3.772, 



ML Relation from Bono et al. (2000 

46.4155 27481.28 

35.6727 19742.83 

28.2406 16894.15 

25.3960 14460.62 

22.3804 12772.75 

19.3886 11263.84 

17.7750 9034.392 

12.6085 5637.467 

10.2940 4371.170 

9.81474 4004.569 

8.37226 3217.688 

8.12498 3032.639 

7.69466 2864.745 

7.23505 2705.430 

6.01283 2105.589 

ML Relation from Chiosi 



5445, 


.00 


18329.74 


4826, 


.00 


5434, 


.01 


14321.20 


4962, 


.65 


5502, 


.35 


12093.99 


4978, 


.72 


5562, 


.28 


10570.70 


5010, 


.85 


5615, 


.34 


9283.257 


5069, 


.43 


5687, 


.98 


8143.460 


5145, 


.20 


5561, 


.39 


7044.743 


5147, 


.04 


5528, 


.15 


4800.840 


5278, 


.11 


5537, 


.42 


3839.001 


5377, 


.80 


5518, 


.39 


3580.396 


5383, 


.61 


5596, 


.34 


2929.354 


5437, 


,17 


5581, 


.53 


2793.441 


5438, 


.96 


5602, 


.37 


2640.439 


5467, 


.54 


5639, 


.03 


2490.271 


5504, 


.98 


5669, 


,57 


1986.586 


5567, 


.91 



40.2561 


23500.98 


5821.99 


11853.38 


4947, 


.11 


35.4374 


19604.60 


5830.18 


9898.124 


4962, 


.26 


28.2629 


14661.29 


5875.98 


7659.896 


5057, 


.96 


25.6378 


12836.04 


5884.01 


6916.299 


5108, 


.02 


19.4314 


9538.226 


5985.21 


5555.179 


5303, 


.60 


17.5523 


8035.700 


5921.07 


4947.967 


5142, 


.58 


14.3143 


5490.345 


5783.01 


3639.964 


5164, 


.11 


11.3659 


3850.971 


5814.49 


2722.473 


5242, 


.47 


10.0011 


3286.030 


5817.00 


2377.950 


5300, 


.52 


9.77393 


3136.139 


5801.77 


2301.751 


5296, 


.73 


8.94637 


2738.003 


5798.30 


2071.106 


5331, 


.19 


8.31297 


2499.416 


5794.49 


1917.696 


5365, 


.34 


8.10751 


2374.713 


5780.35 


1850.048 


5369, 


.40 


8.09994 


2235.292 


5723.18 


1796.586 


5355, 


.41 


7.71463 


2271.393 


5805.49 


1772.270 


5408, 


.34 


ML Relation from 


interpolated Chwsi 


fl989'l 




30.9661 


14958.03 


5611.77 


8486.049 


4963, 


.81 


15.9744 


5697.042 


5665.97 


4209.487 


5079, 


.81 


10.0562 


3197.451 


5670.64 


2560.871 


5283, 


.85 


8.51317 


2723.566 


5716.48 


2201.628 


5382, 


.21 


6.46303 


1852.537 


5684.63 


1593.714 


5507, 


.38 



BC = -0.0153 + 2.122Ar - 0.02 log(g) 
-11.65(Ar)^ 

Note that these functions are also obtained from the BaSeL 
atmosphere database. 

We can compare the colours obtained from the BaSeL 
and SBT atmosphere for our models constructed in this pa- 
per. The results are presented in Figure |K| From this figure 
it can be seen that the colours obtained from both of the 
atmosphere fits agree within 0.05mag. level. The difference 
is even smaller if the micro-turbulence velocity of 5.0kms~^ 
is used. This indicates that the {V — I) colours can be used 
to indicate the temperature. Since the results of our models 
are qualitatively compared to the observations (see next sec- 
tion) and not used to quantitatively derive any theoretical 
PC and/or AC relations, an accuracy of ~ 0.05mag., inde- 
pendent of period, from the atmosphere fit is acceptable and 
would not cause problems for our results. Note that the SBT 
atmosphere are only defined for 5000A' ^ Te// ^ TOOOAT 
and 0.75 ^ log{g) ^ 3.00, few of our models either the Te// 
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Figure 5. Comparison of the colours at maximum, mean and 
minimum light from the models using different atmosphere 
databases. The j;-axis indicates the difference when using the 
BaSeL atmosphere and the SET atmosphere, i.e. BaSeL — SBT. 
The circles, squa r es an d triangl e s are for the models cal culated 
with iBono et alJ i200(Tl) . Ichio^ il98Hh and interpolated IChiosil 
il98gr ML relations, respectively. Open and filled symbols cor- 
respond to the use of micro-turbulence velocity of 1.7kms~^ and 
5.0kms~^, respectively, from the SBT atmosphere database. 

or \og{g) or both are beyond these ranges at certain phases, 
hence no colours can be obtained from the SBT atmosphere 
(for example some points are missing at minimum hght for 
few of the long period models, as shown in Figure]^. Due to 
these reasons, we continue adopt the BaSeL atmosphere fits 
to convert the temperatures and (V — I) colours, after tak- 
ing account of the effective gravity in the fits, as a function 
of phase. 



6 RESULTS FROM THE MODELS 

The effective temperatures for the full amplitude models in 
Table 13 at the corresponding maximum and minimum light 
(or luminosity) are given in Table 0] For the effective tem- 
peratures at mean light, the temperatures for the mean light 
at ascending and descending branch of the light (or luminos- 
ity) curve are not the same (e.g., in Paper II), hence Table|S] 
gives the effective temperature at these phases for our LMC 
models. The layout of Table|21is the same as table 3 from Pa- 
per II. Following Paper II, the locations of the photosphere 
can be identified in the temperature and opacity profiles. 
These are displayed in Figure 1^151 with a log(P) > 1.0, a 
log(P) = 1.0 and a log(P) < 1.0 model, respectively. The 
left and right panels of Figure ISISI are the temperature and 
opacity profiles respectively. The photospheres are marked 
as filled circles in these figures. Finally, the plots of the A, 
the "distance" between the photosphere and the HIF from 
the temperature profiles, as a function of pulsating period 
for the LMC models are presented in Figure|5]with the three 
ML relations used. In Paper II, it is found that the distri- 
bution of A as a function of period is almost independent 



of the adopted ML relation. This is also seen in the LMC 
models as depicted in Figure]^ 

Figure and Figure |5] bear witness to the fact that at 
maximum light, the photosphere lies at the base of the HIF 
for all of the models. Although there is a slight deviation for 
some longer period models, the location of the photosphere 
is close to the HIF within the error bars (which are defined as 
the coarseness of the grid points around the location of the 
HIF). As in Paper II for the Galactic models, the closeness 
of the photosphere to the base of the HIF, for reasonably low 
densities, results in a flat or almost flat PC relation for the 
long period LMC Cepheids. In the case of minimum light, 
even though Figure |^ implies that A{MIN) is nearly con- 
stant across the period range and the photosphere is near the 
base of the HIF, as in the case of maximum light, A{MIN) 
does follow a shallow correlation with period after 10 days. 
Judging from the error bars of A{MIN) and from Figure 
16181 there is tentative evidence that the photosphere is dis- 
engaged from the HIF for log(P) > 1.0 at minimum light. 
Hence the temperatures or the colours at minimum light are 
more dependent on period for log(P) > 1.0 and the global 
properties. 

Theoretical quantities that can be computed from the 
models and compared with data include the pulsation pe- 
riods, the l/-band amplitudes and the Fourier parameters, 
the temperatures and colours at the maximum, mean and 
minimum light. These are the PC plots, the AC plots, the 
period-temperature plots and the Fourier parameters plots 
portrayed in Figures I10I13I The temperatures in Table 
|S] after conversion to the {V — I) colours as mentioned in 
previous section, are superimposed along with the observed 
LMC PC relations as plotted in Figure 1101 Similarly, Fig- 
ure llll graphs the same quantities but on the log(r)-log(P) 
plane with the observed {V — I) colours converted to tem- 
peratures using the prescriptions given in Section 5. The 
theoretical bolometric light curves are converted to the V- 
band light curves with the bolometric corrections obtained 
from the BaSeL database mentioned previously. From the 
theoretical V-band light curves, the amplitudes can be esti- 
mated and these are displayed in Figure WH along with the 
colours from models to compare with the empirical AC re- 
lations. The Fourier parameters of the theoretical l/-band 
light curves can also be obtained with (n — 6) Fourier ex- 
pansion. These Fourier parameters are compared with the 
observational data in Figure IT!T1 Several features are noticed 
from Figure colli 

(i) The general trends of the models qualitatively match 
the observational data. There are greater discrepancies be- 
tween the data and short period models, particularly in 
matching the observed light curve amplitude s. 

(ii) The models with the ML relation from lChiosil Jl989l) . 
with lower L/M ratio, do better in matching the observa- 
tions. These models also tend to lie near the envelopes of the 
PC, AC, log(r)-log(P) and ^„(F)-log(P) relations defined 
by the observational data. 

(iii) The slopes of the period-colour (or period- 
temperature) relations at maximum and minimum light 
from the models roughly match the observational data, i.e., 

* The slopes of the A-log(P) relation may not be correlated with 
the slopes of the PC relation 
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Table 5. Temperatures at mean light from full-amplitude non-linear model calculations. See Paper II for the meanings of < L >, L 
Tmean and T . The periods, luminosity and temperature are in days, Lq and K, respectively. 
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<L> 


Lmean {dSCj 
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ML Relation from Bono et al 
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the theoretical PC(max) relation is approximately flat and 
there is a relation at minimum light. 

(iv) The temp eratures from the models with the 
iBono et alJ (I2OO0I) ML relat i on is c ooler (hence redder) than 
the models with the lChiosil I^S^ ML relation and the ob- 
served data at maximum light. In contrast, the tempera- 
tures (or the colours) at minimum hght from the models 
with these two ML relations are consistent with each other 
and are located near the blue edge of the observed data. 

(v) The means at the descending branches are in better 
agreement with the observed data than the means at the 
ascending branches. This is because the observed means, 
{V — j'l'P^'^^^^^ are obtained mostly from the descending 
branches. Though previous researchers have noted that tem- 
peratures on the ascending and descending branches are not 
the same at mean light (as Cepheids exhibit loops in CMD), 
what is new here is the way the nature of the HIF changes 
during the pulsation. 

(vi) The behaviors of the models from the interpolated 



'Chiosi ("1989^ ML are closer to the models from lBono et alJ 
( 2000) ML relation because their slopes are very similar. 

(vii) The amplitudes of the theoretical light curves (in 
both of the bolometric and V^-band hght curves) are smaller 
than the obser vations at given p eriod, especially for the 
models with the lBono et alJ ll200(]|) ML relation. These can 
be seen from the AC relations as given in Figure IT^ and the 
left panels of Figure [1^1 

Overall, some agreements and disagreements are found be- 
tween the theoretical quantities and the observational data. 
It is also found out that there are some problems associated 
with the pulsation codes when the LMC models are con- 
structed: these include the smaller amplitude of the model 
light curves and the cooler temperatur es at the maximum 
light (especially with iBono et alJ I2OOOI ML relation). Note 
that from equation (1), cooler temperatures at maximum 
light imply that the amplitudes will be lower at given period. 
Varying other parameters in the pulsation codes, including 
the a parameters, does not improve the situation, though 
perhaps a more detailed and systematic study of the depen- 
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Figure 6. The temperature (left panel) and the opacity (right panel) profiles for a long period LMC model. The dotted, solid and dashed 
curves are for the profiles at maximum, mean and minimum light, respectively. The filled circles are the location of the photosphere at 
T = 2/3 for each phases. The mean light profiles at the ascending and descending branch are the solid curves that lie close to the profiles 
at maximum light (dotted curves) and minimum light (dashed curves), respectively. 
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Figure 7. Same as Figure |H] but for a 10-days period LMC model. 



dence of LMC Cepheid pulsation models on the a param- 
eters could resolve this situation. However, we believe that 
the qualitative nature of the photosphere-HIF interactions 
as given in Figure]^ will still hold even in models which fare 
better in mimicking observed amplitudes. This is in part be- 
cause Figure|U]suggest that the behaviors of A as a function 
of period are nearly, though not completely, independent of 
amplitudes, as the models with Chiosi (1989) ML relation 
have higher amplitudes (although st ill smaller than th e ob- 
servations) than the models with the lBono et alJ i200Ci) ML 
relation. However, better codes that fix these problems or 
the 3-D convection codes are needed in the future studies. 



6.1 Comparison with the Galactic Models 

The temperature profiles from the Galactic models given in 
Paper II and the LMC models are compared in Figure ITU at 
maximum and minimum light. The upper panels of Figure 
1141 suggest that at maximum light, the photosphere is not 
far from the base of the HIF in both of the Galactic and the 
LMC models. In contrast, the photosphere is further away 
from the HIF in the Galactic models than the LMC mod- 
els at minimum light. The HIF is located further out in the 
mass distribution for the Galactic models. The plots of the 
A-log(P) relation from the Galactic and LMC models at 
maximum and minimum light are also compared in Figure 
1151 It can be seen from the figure that at maximum light, the 
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Figure 8. Same as Figurel^l but for a short period LMC model. 
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Figure 9. The p lots of A as function of lo g(P). The open squares, solid square s and open circles are the models calculated with 
iBono et alj 120001 ML relation, IChiosil ilQSgtl ML relation and interpolated IChiosil ll98Sll ML relation, respectively. The dashed lines 
represent (roughly) the outer boundary of the HIF. 



behavior of both Galactic and LMC models is similar, where 
the photosphere is near the base of the HIF. At minimum 
light, the long period models show that the photosphere is 
disengaged from the HIF, while the behavior of the short 
period models is different between the Galactic and LMC 
models. The photosphere of the short period LMC models 
seems to be located closer to the HIF at minimum light, but 
it is not the case for the short period Galactic models. This 
could lead to shallower slopes of the PC(min) relation seen 
in the LMC Cepheids as compared to the Galactic counter- 
parts. In terms of the HIF-photosphere interaction, there is 
some tentative evidence from the models that the LMC long 
period Cepheids behave like the Galactic Cepheids, while the 



short period LMC Cepheids behave like the RR Lyrae stars 
at minimum light. 

Figure IHH graphs the density (defined as 1/V, where V 
is the specific volume) at the photosphere as a function of the 
period of the model at minimum, maximum and ascending 
and descending mean light. Galactic models generally tend 
to have the lowest density and, in particular, have signifi- 
cantly lower densities at minimum light than the LMC mod- 
els. We note that the Galactic models always have a pho- 
tospheric density lower than about 10~*g/cm^ whereas the 
photospheric density for the LMC models only falls below 
this figure after a period of 10 days. At maximum light, all 
long period models have a low photospheric density. What 
we get from this figure is that it provides some evidence 
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Figure 10. The LMC PC relations with the results from the models. The smal l crosses are the observed data points as given in Figure^ 
The open squares, solid squares and open circles are the models calculated with lBono et al ] HOOO), Chiosi (1989) and interpolated Ch iosil 
jl989() ML relations, respectively. The temperatures of the models are converted to the {V — I) colour using the BaSeL database. Left 
panel: PC relations at maximum and minimum light. Right panel: PC relations at mean light for both of the ascending and descending 
means. The mean colours for the observed data are the {V — /)P'^™'^^" as given in Figure 1^ 
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Figure 11. The plots of log(T)-log(P) relations for the LMC data and models. The symbols are the same as in Figure lTol The conversion 
of the (V — I) colours to the temperature are done using the equations given in Bcaulieu et al. |2001j). Left panel: log(T)-log(P) relations 
at maximum and minimum light. Right panel: log(T)-log(P) relations at mean light for both of the ascending and descending means. 



that there is a difference in photospheric density between 
the LMC and Galactic models. Moreover, this difference ap- 
pears to be consistent with what is required by our theoreti- 
cal scenario: short period LMC models have a higher photo- 
spheric density than their Galactic counterparts. However, 
for a discussion of some caveats, see Section 8. 



7 APPLICATIONS 

We now discuss two important applications of the 
photosphere-HIF interaction: reddening corrections and the 
explanati on of t he ob served non-linear LMC PL (and PC) 
relations. ICodd 1^43) original interest in the spectral prop- 
erties of Cepheids at maximum light was to estimate red- 
dening. SKM used th is to correct a number of reddening 
for Galactic Cepheids. iFernid il994ft used equation (1) and 
the theoretical explanation provided in SKM to derive a re- 
lation linking the colour excess to the colour at maximum 
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Figure 12. The LMC AC relations with the results from the models. The small crosses arc the observed data points as given in Figure 
m Th e open squares, solid squares and open circles are the models calculated with Bono et al. (2000), Chiosi (1989) and interpolated 
IChiosil Il98df) ML relations, respectively. The bolometric light curves from models are converted to V-band light curves with the BC 
obtained from the BaSeL database. Left panel: AC relations at maximum and minimum light. Right panel: AC relations at mean light 
for both of the ascending and descending means. The mean colours for the observed data are the {V — /)P'^™'^^" as given in Figure HI 
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Figure 13. The Fourier parameters as function of periods for the LM C obse rvational data (small crosses) and the models, where open 
squares, solid squares and open circles are the models calculated with iBono et al. (2000), Chiosi 1 1989) and interpolated Chiosi 1 19i^ 
ML relations, respectively. Upper Left: Plot of A\ vs log(P) in V-band. Lower Left: Plot of A2 vs log(P) in V-band. Right: Plot of if>2i 
vs log(P), where 4>2i = 4>2 - 2<^i <Simon fc Ledll98lt iNgeow et al.ll2003t) . 



light, the V-band amphtude an d the p e riod. Such a relation 
is predicted from equation (l). iFeriiig lll994l) estimates the 
error with this method to be comparable to other multi- 
colour methods. 

A more interesting application of the HIF-photosphere 
interaction is to explain the recent detected n o n-line ar LMC 
P L relation as presente d in iTammann et al] ||2002|) . Paper 
I, ISandaee et al.l t004) and et al.l j2005lP Paper I 

used the F-test to provide strong statistical evidence that 
the optical Cepheid PL relation at mean light in the LMC 



is non -linear around a period close to 10 davs. lNgeow et alJ 
f2005) used the MACHO and 2MASS datasets together with 
additional long period Cepheids from the literature to fur- 
ther support the existence of non-linear LMC PL relation in 
the optical and near infra-red wave-bands. In contrast, cur- 
rent data in dicate that the Galactic PL relation is linear a t 
mean light jTammann et al.ll2003l: iNgeow fc Kanbuill2004l) . 

Non-linearity of the LMC PL relations can be tested 
using the F-test with the data given in Section 3. The 
empirical results of the fitted LMC PL relations at max- 
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Figure 14. Comparisons of the temperature profiles between the Galactic and LMC models at maximum (upper panel) and minimum 
(lower panel) light. The solid and dashed curves represent the temperature profiles from the Galactic and the LMC models, respectively. 
The left and right panels are for the short and long period models, respectively, with the periods labeled at the upper left corners. The 
filled (for the Galactic models) and open (for the LMC models) circles are the locations of the photosphere in these models. 
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Figure 15. Comparisons of the A-log(P) plots between the Galactic (left panel) and the LMC (right panel) models at maximum (top 
panel) and minimum (bottom panel) light. 



imum, mean and minimum light using the updated data 
are presented in Table |S| The plots of the PL relations at 
maximum/minimum light and at mean light are shown in 
Figure 1171 & 1181 respectively. The F-test results for these 
PL relations are: Fv (max, mean, min) = {2.03, 8.22, 16.1}, 
and (max, mean, min) = {0.28, 7.37, 17.9}. The large F- 
values for both V- and 7-band PL relations at mean and 
minimum light strongly indicate that the PL relations at 
these two phases are not linear, and the data is better de- 
scribed with the broken (i.e, two regressions) PL relation. 
However, the small F- values at maximum light, with corre- 
sponding p-values of 0.13 and 0.76 for the V- and 7-band 



PL(max) relations respectively, show that the null hypoth- 
esis of the F-test cannot be rejected (a value of p < 0.05 
and/or f > 3 is required for doing this). Hence there is no 
observed break seen in the PL(max) relation and the data 
is consistent with single line regression. Note that the same 
slopes of the PL(max) relations for long period Cepheids in 
both bands are consistent of the finding that the PC(majc) 
relation is flat for these Cepheids. 

Our tentative theoretical explanation for the non-linear 
nature of the LMC PL relations across a p eriod of 10 days 
repl ies on the HIF-photosphere interaction. ISandagj 1^53) 
and iMadore fc FreedmanI (Il99 j) have established the con- 
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Figure 16. Log of the photospheric density plotted against period for tlie pulsation models. The symbols are: open squares = LMC 
models with Bono ML, filled squares = LMC models with Chiosi ML; open circles = LMC models with interpolated Chiosi ML; crosses 
= Galactic models with Chiosi ML. 



Table 6. The period-luminosity relation in the form of m^^jj = a(v,/) ^og{P) + ^'nd '^(v,l) dispersion of the relation. 



Phase ay by ay clj bj aj 















All, N = 641 
















Maximum 


-2.792 


± 


0.044 


16.706 


± 


0.031 


0.260 


-2.966 


± 


0.029 


16.380 


± 


0.020 


0.170 


Mean 


-2.735 


± 


0.035 


17.088 


± 


0.025 


0.208 


-2.963 


± 


0.024 


16.609 


± 


0.017 


0.141 


Minimum 


-2.546 


± 


0.034 


17.270 


± 


0.024 


0.204 


-2.837 


± 


0.023 


16.738 


± 


0.017 


0.140 














Long 


Af = 63 
















Maximum 


-3.014 


± 


0.180 


16.959 


± 


0.225 


0.257 


-3.040 


± 


0.122 


16.469 


± 


0.152 


0.174 


Mean 


-2.698 


± 


0.160 


17.099 


± 


0.199 


0.228 


-2.982 


± 


0.111 


16.671 


± 


0.139 


0.158 


Minimum 


-2.315 


± 


0.171 


17.052 


± 


0.213 


0.243 


-2.617 


± 


0.122 


16.512 


± 


0.153 


0.174 














Short, 


N = 578 
















Maximum 


-2.680 


± 


0.077 


16.640 


± 


0.048 


0.260 


-2.946 


± 


0.050 


16.368 


± 


0.031 


0.169 


Mean 


-2.937 


± 


0.060 


17.205 


± 


0.038 


0.203 


-3.090 


± 


0.041 


16.683 


± 


0.026 


0.138 


Minimum 


-2.818 


± 


0.057 


17.429 


± 


0.036 


0.194 


-3.030 


± 


0.039 


16.852 


± 


0.024 


0.131 



nection between the PC and PL relations; both these rela- 
tions arise from the more general PLC relation. These rela- 
tions refer to quantities evaluated at mean light. The exis- 
tence of such a connection relies on the period-mean density 
theorem, the instability strip and the Stefan-Boltzmann law. 
If we assume the Stefan-Boltzmann law can be applied at 
every phase, then it is straightforward to show that a PLC 
relation (though possibly with different coefficients) exists 
at every phase point. Thus the standard PLC relation and 
indeed the PC and PL relation expresses at mean light are 
just the averages of the same relations at different phases 
points. Consequently one way to understand the behavior 
of PLC/PL/PC relations at mean light is to understand 
their behavior at different phase points. What we try to 
do in this paper is point out some evidence from our models 
that shows how the changing behavior of the PC relations 
at different phases can, in principle, arise from a consider- 
ation of the photosphere-HIF interaction at these phases. 
Since the mean hght PC and PL relation are the average of 
those at all phases, these properties can affect the PC and. 



as a consequence, the PL relation (via the PLC relation). 
In fact, the new data with superb phase resolution from 
such micro-lensing projects such as OGLE and MACHO de- 
mands a multiphase analysis. This approach can potentially 
lead to a deeper understanding of the pulsation and evo- 
lution of Cepheid variables. For example. iNgeow fc Kanbuil 
(2006) looked at PC relations in the Galaxy and LMC as 
a function of phase. They found that short and long period 
LMC Cepheids have a shallower and steeper slope at most 
pulsation phases than Galactic Cepheids respectively. 



8 CONCLUSION AND DISCUSSION 

In this paper, we have confronted updated PC and AC re- 
lations at maximum, mean and minimum light for LMC 
Cepheids observed by the OGLE team, and additional 
Cepheids from the literature, with theoretical, full amplitude 
pulsation models of LMC Cepheids. The observed PC and 
AC relations provide compelling evidence of a non-linearity 
or break at a period of 10 days. We also constructed theoreti- 
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Figure 17. The LMC PL relation at maximum (left panel) and minimum (right panel) light. The upper and lower panels are for the 
y-band and the /-band PL relations, respectively. The open and filled circles are for short and long period Cepheids, respectively. The 
dashed lines are the fitted PL relations for the long period Cepheids. While the solid lines are the fitted PL relations for the short period 
Cepheids, and extended to the longer period range to compare with long period PL relations. 
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Figure 18. Same as Figure ITtI but for the LMC PL relations at 
mean light. 

cal Cepheid pulsation models appropriate for th e LMC using 
the Florida pulsation codes (lYecko et aLlligO^ to study the 
HIF-photosphere interaction. 

The empirical results pres ented in this p a per, a s well 
as in other pape rs such as ISandaee et alJ (|20^ and 
iNgeow et alJ (l200,'j ) . provide strong empirical evidence that 
the PC and PL relations for the LMC Cepheids are non- 
linear, in the sense described in previous sections. Issues 
such as extinction and a lack of long period Cepheids that 
may cause the non-linear LMC PL and PC relations have 
been a ddres sed and argu ed against in Paper L iSandage et alJ 
l)2004) and iNgeow et a l. (2005), and will not be repeated 
here. Other arguments against the non-linear LMC PL re- 



lation include the results presented in lPersson et alJ i2004) . 
as the authors found no evidence for a non- linear PL re- 
lation in the LMC at JiJ J^-bands. However. [Ngeow et alJ 
(2005) treated the data of iPersson et al.| i2004l) extensively 
and found, in a st a tistica lly rigorous way, that the reason 
whv lPersson et all (|20o3) found linear JHK PL relations, 
is due to the sm all number of s hort p eriod Cepheids 18) 
in their sample. iNgeow et alj (l2005l) also reduce the num- 
ber of OGLE/MACHO LMC Cepheids and show how the 
F-test can produce a non-significant result when the num- 
ber of short/long period Cepheids become small. Instead, 
using the 2MAS S data that a r e cros s-correlated with MA- 
CHO Cepheids, iNgeow et all have found that the 
LMC J/f-band PL relations are no n-linear^ and the K- haxiA 
PL relation starts to become linear. iNgeow et alj i2005l) also 
discussed why this is the case. Another argument against 
the non-linear PL rela tion is that th e PL r elation should be 
universal, as found in ICieren et alj ll2005l) . We argue that 
their results are based on a handful of Cepheids (~ 15) 
and on short periods Cepheids in a cluster whose member- 
ship to the LMC is in question. Their shallower Galactic 
PL relation based on the revised infra-red surface bright- 
ness method also contradicts the steeper Galactic PL rela- 
tion based on ind ependent methods from open cluster main- 
sequcncc fitting ( Tam mann et al] l200.?l : iNgeow fc Kanbuil 
i2004; Sandagc ct al. 20oJ! 

It is worthwhile to point out that our sample selection 
does not affect the detection on non-linear LMC PL rela- 
tion at mean light. Since the mean magnitudes of a Cepheid 
light curve is less affected by our constrains on selecting 
the Cepheids with good light curves, we can use the pub- 
lished (reddening corrected) mean magnitudes to test the 
non-linear LMC PL relation. The anonymous referee kindly 

^ Non-linear PL relations a re relatively easier t o see in V-band 
than in J-band, as shown in lNgeow et alJ i2005l) . 
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provided a large sample of LMC Cepheids that combined 
the published me an F -band magnitudes f rom t he OGLE, 
ISebo et alJ ^o3) and ICaldwell fc LanevI I^1991^ datasets. 
There are a total of 115 long period Cepheids in this sam- 
ple and the F-test still return a significa nt detection of the 
non-linear LMC PL relation. The OGLE-f jSebo et al.l (l2002l) 
combined data also g ive very similar res ults. Similar tests 
have also been done in lNgeow et al. ( 2005") by using the MA- 
ClfO data alone and the MACH04 Scbo c t al. (JOQgj com- 
bined data. The non-linear LMC PL relation is still present 
from the _F-test results on these two datasets. Therefore we 
believe our sample selection does not affect the detection of 
the non-linear LMC PL relation. The detection of non-linear 
LMC PL relation from totally independent OGLE and MA- 
ClfO data, using totally independent reddening estimates, 
suggested that this non-linearity is real and our paper is the 
first attempt to theoretically explain this non-linearity in 
terms of the HIF-photosphere interaction. 

Due to small number of LMC models, it is impossible 
to derive the theoretical PC and AC relations with a small 
error on the slope and compare directly to the empirical 
relations. However, these LMC models can be qualitatively 
compared to the observations by converting some physical 
quantities to the observable quantities and vice versa, such 
as the temperature-colour conversion. Hence we compared 
our model light curves to the observations in terms of the- 
oretical PC and AC relations at the phases of maximum, 
mean and minimum light and also in terms of the Fourier 
parameters from theoretical light curves with observations. 
The theoretical quantities from the models generally agree 
with the observations, but it was found out that these mod- 
els tend to have smaller amplitudes and (hence) the tem- 
perature is cooler at maximum light than the real Cepheids. 
Though our models have some drawbacks in this compari- 
son, our main interest is in comparing the interaction of the 
photosphere and HIF as a function of phase with similar re- 
sults presented in Paper H for Galactic Cepheid pulsation 
models. The aim is not to compare our models rigorously 
with observations but rather to study models which match 
observations reasonably well in the context of the theoreti- 
cal framework described in previous sections and in Paper I 
& H. Nevertheless we argued that the qualitative nature of 
the photosphere-HIF interaction is not seriously affected by 
these problems. 

Our postulate is that at certain phases, this interac- 
tion can affect the PC relation due to the properties of the 
Saha ionization equation: specifically for reasonably low den- 
sities in Cepheid envelopes, hydrogen ionizes at a temper- 
ature that is almost independent of period. Consequently, 
when the photosphere is located at the base of the HIF, the 
photospheric temperature and hence the colour is almost 
independent of period. However, when this engagement oc- 
curs, but the density is greater, then the temperature at 
which hydrogen ionizes again becomes sensitive to global 
surroundings and hence on period. When the photosphere 
is not engaged with the HIF in this way, its temperature is 
again dependent on period and global stellar parameters. 

For Galactic Cepheids, this HIF-photosphere interac- 
tion occurs mainly at maximum light for Cepheids with 
logP > 0.8 (Paper II). At minimum light, there is a strong 
correlation between the HIF-photosphere distance and pe- 
riod leading to a definite AC relation at minimum light for 



Galactic Cepheids (SKM, Paper I & II). In this paper, we 
have found tentative evidence that, for short period LMC 
models which match observations in the period-color plane, 
the HIF-photosphere interaction occurs at most phases but 
at densities which are too high to produce a flat PC re- 
lation. Why would these short period LMC Cepheids be 
different in this regard to short period Galactic Cepheids? 
One possibility could be that this is partly because these 
LMC Cepheids are hotte r than their Ga l actic counterparts 
jKanbur fc Ngeowl l2004l: ISandaee et ai]l2004) . The HIF- 
photosphere are disengaged for most of the pulsation cycle 
for long period LMC Cepheids. This happens because as the 
period increases, so does the L/A'I ratio which pushes the 
HIF further inside the mass distribution. When the HIF- 
photosphere are disengaged in this way, the photospheric 
temperature is more dependent on density and hence on pe- 
riod. The change is sudden because the HIF-photosphere 
are either engaged or they are not. This can lead to a sud- 
den change i n the PC relation at 10 days as shown by the 
observations lE ^mmam^^l .12002 : Kanbu r fc Ngeowl2004l: 
ISandage et al.ll2004HNgeow et al.ll2005,) . However, at maxi- 
mum light the HIF-photosphere are engaged at low densities 
for long period LMC Cepheids leading to the observed flat 
PC relation for these stars. Taken together with equation 
(1), this theoretical scenario is consistent with the observed 
PC- AC behavior described in Paper I and in this study. The 
anonymous referee has noted that these suggestions about 
photospheric density can be tested by spectroscopic means. 

We now enumerate some caveats to our argument that 
could be addressed in future papers. 

(i) Since the SMC PC relation at mean light is linear 
(e.g.. Paper I), how do SMC (i.e., metal-poor) models fit 
into the theoretical scenario outlined in this paper and Pa- 
per II, if at all? This is a difficult question and its full answer 
is beyond the scope of this paper. However, as the metal- 
licity decreases, we do note that the SMC has a different 
ML relation to the LMC and Galaxy and so does the tem- 
peratures associated with the instability strip. These will 
change the re l ative location of the HI F and photosphere 
jKanbuilll995l: iKanbur fc Phillip j 119961) and possibly alter 
the phase at which they interact. Further the amplitudes 
for SMC Cepheids are sma ller due to the lower metallicity 
iPaczvhski fc Pindoi<l20o3) . This will also affect the HIF- 
photosphere interaction. One difference which can be con- 
sistent with this is the fact that the PC relation at maximum 
light in the SMC is not flat (see Paper I) but it is the case 
for the Galaxy and LMC PC relations. This indicates that 
at maximum light, there is less interaction between the HIF 
and photosphere at low densities. This leads to an observed 
linear PC relation at mean light for the SMC Cepheids. 
These will be investigated further in a future paper in this 
series. 

(ii) Could the well-known Hertzsprung progression play 
any part in causing the observed changes in the Galactic 
and LMC PC relations? 

(iii) It may also be that higher order overtones becoming 
unstable or stable, though with the fundamental mode still 
being dominant, may also have an impact on the PC relation 
in some as yet unknown way (Paper II). 

(iv) The behavior of short period LMC Cepheids still 
needs to be understood, for example, what causes the dif- 
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ference between the bottom left panels of Figures El and 
I15f That is, why is it that for short period Galactic/LMC 
Cepheids, the HIF-photosphere are disengaged/engaged? 
Our experience suggests that constructing short period full 
amplitude fundamental mode Cepheids requires more care 
than the long period case because the first overtone has a 
non-negligible growth rate. Because of this we feel a thor- 
ough study of these short period Cepheids merits a separate 
paper. 

(v) Would more advanced pulsation codes which, for ex- 
ample, can match the observed amplitudes and which con- 
tain a more accurate model of time dependent turbulent 
convection, yield similar results, especially for Figure 1151 :' 
Could such codes fare better in modeling short period LMC 
Cepheids? 
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